1 



o 
o 

> 

o 

00 



o 

CD 



O 

o 



in 



o 



X 



Comment on "Phase Transitions in Systems 
of Self-Propelled Agents and Related Network 
Models" 

In a recent letter [1| , Aldana et al. study order-disorder 
phase transitions in random network models and show 
that the nature of these transitions may change with the 
way noise is implemented in the dynamics. Arguing that 
these networks are limiting cases of simple models of in- 
teracting self-propelled agents of the type of the Vicsek 
model (VM) [3| , they claim that the conclusions reached 
for the networks may carry over to the transitions to col- 
lective motion of the VM-like systems. They suggest in 
particular that in the case of "angular" noise (i.e. as in 
the original VM or in their Eq. (1)) the transition 
to collective motion is continuous^ in contradiction with 
some of the conclusions of Q. While we agree with the 
analysis of the network models, we argue here that it has 
no bearing on VM-like systems. We show in particular 
that the transition to collective motion, for angular noise, 
remains discontinuous for any finite microscopic velocity 
V and finite density p, however large, confirming 

In is was shown that the transition in the origi- 
nal VM appears continuous only when the linear system 
size L is smaller than some crossover size L*{p,v) and 
thus is discontinuous in the thermodynamic limit (see 
Fig. 2a there). These results were obtained for values of 
p and V of order unity. The same scenario occurs in the 
large-w or large p limit of interest here. Figure [1^ here 
displays results akin to Fig. 1 of [H except that the VM 
order parameter curves were obtained in systems such 
that L > L*{p,v), whereas in [H L ~ 100 < L*{p,v) (see 
below for our estimates of L*). Clearly, the transition is 
discontinuous, as demonstrated by the minimum shown 
by the Binder cumulant (Fig. [T]d) and by the bimodal- 
ity of the order parameter distribution in the transition 
region (not shown). 

The key difference between the network models and 
interacting self-propelled agents is indeed that the latter 
move, inducing a local coupling between order and den- 
sity, which is well-known to be crucial for understanding 
collective properties of active particles In the VM, 
this coupling gives rise to strong density and order vari- 
ations on lengthscales of the order of L* . While the net- 
work models in [l[ capture the long-range interactions 
due to large velocities, they obviously cannot account for 
any coupling between density and order. The network 
models only represent VM-like models of size L < L* . 

The crossover scale L* is difficult to estimate with high 
accuracy, but our data indicate that L* increases roughly 
linearly with v (for p = 2, L* ^ 150 ± 25, 175 ± 25, 250 ± 
50, 550 ± 50 for v = 5, 10, 20, 40 respectively). Thus, we 
expect L* to be finite at any finite v. As a consequence, 
the transition is always discontinuous in the thermody- 
namic limit, although its asymptotic behavior is harder 
to observe as v is taken larger and larger. When v is 
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FIG. 1: (Color online) Order-disorder transition for the Vic- 
sek model at different i; values (solid lines) and the network 
of fixed vectors of [l| (dashed line), both with angular noise. 

(a) Time-averaged order parameter (?/)) vs noise strength r;. 

(b) same as (a) for the Binder cumulant G = 1 — {ip'^)'^ / {ip'^) . 
For the VM, p — 2, and v — 5 (circles), i; = 10 (squares), 
V — 20 (triangles) with L = 200, 250 and 300 respectively. 
For the network, the number of nodes is TV = 4 x 10^ and the 
average connectivity is set to K = 2n (i.e. each node inter- 
acts with 6 or 7 neighbors, with suitable probabilities), which 
corresponds to p = 2 in the VM. 



taken to infinity first, the transition is continuous for all 
finite "size" , but then the notion of distance in physical 
space is abolished. 

Summary: the transition to collective motion in VM- 
like systems with angular noise remains discontinuous for 
large v values. Thus, the networks studied in [l[ at best 
constitute a singular v oo limit of these systems Q . 
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In the large p limit, wo obtain similar rosnlts: for L > L* increases with p (results obtained at v = 0.5, not shown), 

the transition reveals its discontirmous character and L* 



